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SrEr 2 C >4 is a geometrically frustrated magnet which demonstrates rather unusual properties at 
low temperatures including a coexistence of long- and short-range magnetic order, characterized 
by two different propagation vectors. In the present work, the effects of crystal fields (CF) in this 
compound containing four magnetically inequivalent erbium sublattices are investigated experimen¬ 
tally and theoretically. We combine the measurements of the CF levels of the Er 3+ ions made on 
a powder sample of SrEr 2 04 using neutron spectroscopy with site-selective optical and electron 
paramagnetic resonance measurements performed on single crystal samples of the lightly Er-doped 
nonmagnetic analogue, SrY 2 0 4 . Two sets of CF parameters corresponding to the Er 3+ ions at the 
crystallographically inequivalent lattice sites are derived which fit all the available experimental data 
well, including the magnetization and dc susceptibility data for both lightly doped and concentrated 
samples. 

PACS numbers: 71.70.Ch, 75.10.Dg, 76.30.Kg, 78.70.Nx 


I. INTRODUCTION 

The family of rare-earth (RE) strontium oxides, 
Sri? 2 04, has recently been identified as a useful addition 
to a small group of compounds where the combination 
of geometrical frustration, magnetic low-dimensionality 
and single-ion physics results in the stabilisation of highly 
unusual ground states at temperatures much lower than 
those expected from the strength of the magnetic inter¬ 
actions [lj, | 2 | . A wide variety of unconventional magnetic 
properties observed in these compounds include a coex- 
istance of long-range antiferromagnetic and short-range 
incommensurate order in SrEr 2 0 4 [H, Q, two types of 
short-range order in SrHo2 04 SB , noncollinear order 
in SrYb20 4 Q, incommensurate magnetic structure in 
SrTb 2 04 @ and the absence of the longer-range mag¬ 
netic correlations in SrDy 2 0 4 down to the lowest ex¬ 
perimentally available temperatures 0 - To elucidate 
and model specific magnetic properties of these oxides, 
which can be described as geometrically frustrated multi- 
sublattice magnets with a quasi-one-dimensional space 
structure, it is necessary to establish the spectral prop¬ 
erties of the magnetic subsystem formed by the RE ions. 
However, only sparse information on low-energy elec¬ 
tronic excitations in Sri? 2 04 compounds has been ob¬ 
tained from studies of inelastic neutron scattering in 
SrHo 2 0 4 mm and SrDy 2 04 0 . The present work is 
aimed at the determination of sets of parameters used in 
the microscopic models of spectral and magnetic prop- 
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erties of a single RE ion in a dielectric crystal. With 
these parameters, it is possible to carry out a physi¬ 
cally justified analysis of the different quantum effects 
caused by the interactions between the paramagnetic RE 
ions. To characterize the single ion properties, we start 
from experimental and theoretical studies of a strongly 
diluted RE subsystem in the isostructural diamagnetic 
compound, SrY 2 04- Studies of SrY 2 C>4 single-crystals 
doped with the Er 3+ ions are well suited for this pur¬ 
pose due to the rich spectrum of electronic excitations 
of the Er 3 + ion lying in the spectral range accessible to 
site-selective optical measurements. 

We have performed a series of spectroscopic studies 
of the Er 3+ ions in the SrY 2 04:Er 3+ ( 0.5 at.%) sam¬ 
ple (measurements of the electron paramagnetic reso¬ 
nance, EPR, and optical spectra) and in the concentrated 
SrEr 2 04 crystals (measurements of the inelastic neutron 
scattering at different temperatures) as well as the mag- 
netometry studies (measurements of the magnetization 
versus external magnetic fields applied along the prin¬ 
cipal crystallographic axes at different temperatures). 
Analysis of the results of these measurements in the 
framework of a semi-phenomenological crystal field (CF) 
theory [0 and the derived four-particle self-consistent 
model of magnetic interactions have allowed us to deter¬ 
mine the sets of CF parameters and exchange coupling 
constants related to the specific positions of the Er 3+ ions 
in the crystal lattice. From comparisons of the simulated 
spectral envelopes of the inelastic neutron scattering, the 
temperature dependencies of the dc susceptibility tensor 
components and the field dependencies of the magnetiza¬ 
tion in SrEr 2 04 with the experimental data, it is evident 
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that the spectral parameters of the impurity Er 3+ ions 
in SrY 2 04 can be used successfully to reproduce spec¬ 
tral and magnetic properties of the erbium subsystem in 
SrEr 2 C> 4 . 

The paper is organized as follows: the experimental 
methods employed in the present work and the crystal 
structure of the samples are described in Sec. EH The 
results of measurements are presented in Sec. EIH Sec.EY ] 
contains an analysis of the data obtained which involves 
simulations of the crystal fields affecting the Er 3+ ions 
at different crystallographic sites and modeling of the 
magnetic properties of the concentrated SrEr 2 C >4 crystals 
in the paramagnetic phase. The paper ends with a brief 
summary of the results. 



II. EXPERIMENTAL METHODS AND 
CRYSTAL STRUCTURE 

A single crystal of SrY 2 C >4 doped with Er 3+ ions (0.5 
at.%) was grown by the optical floating zone technique. 
SrC 03 (Alfa Aesar, 99.99%), Y 2 C >3 (Alfa Aesar, 99.99%) 
and Er 2 03 (Alfa Aesar, 99.9%) oxides were used as start¬ 
ing materials. The thoroughly ground and mixed stoi¬ 
chiometric composition was initially sintered at 1050 °C 
for 8 hours. The synthesized powder was examined with 
powder X-ray diffraction (Bruker D8 ADVANCE, Cu 
Ka) and found to be a phase-pure material with lattice 
parameters and space group symmetry corresponding to 
the structure of SrY 2 C> 4 . A cylindrical rod was formed 
from the powder using a hydrostatic press, and this rod 
was used as the ingot for the crystal growth. The single 
crystal was grown in an air flow of 0.5 L/min at ambi¬ 
ent pressure using an optical floating zone furnace FZ- 
T-4000-H-VII-VPO-PC (Crystal Systems Corp., Japan) 
equipped with four 1 kW halogen lamps at a rate of 
3 mm/h. The feed and seed rods were counterrotated 
at a rate of 15 rpm in order to obtain a homogeneous 
molten zone. The grown crystal was pinkish in color and 
was of high optical quality. 

The crystal structure of Sri? 2 C >4 is orthorhombic with 
the space group Pnam. The unit cell has 28 ions in 
total, split into seven groups each containing four crys- 
tallographically equivalent A-ions (A = Sr, Rl, R2, 01, 
02, 03 and 04) at the Wyckoff 4c positions. The ions 
coordinates (in units of lattice constants a, b, and c), 
ri, a = (x\, y\, 0.25), r 2jA = (0.5 - x x ,y\ - 0.5,0.75), 
Us,A = {—xx, —yx, 0.75), r 4jA = (x\ - 0.5, 0.5 - y x , 0.25) 
are determined by structure parameters x\ and y\. The 
lattice constants and the structure parameters of SrY 2 04 
and SrEr 2 04 were reported in Refs. [ll. fill. Il5j|. In par¬ 
ticular, a = 10.0166 A, b = 11.8541 A, and c = 3.3848 A, 
a-'Ri = 0.4227, j/ R i = 0.1101, and a; R2 = 0.4232, y R2 = 
0.6118 in SrEr 2 04 at 80 K [l5j. The impurity Er 3+ ions 
in SrY 2 C >4 substitute for Y 3+ ions at the Rl and R2 sites 
and, as in SrEr 2 C> 4 , have six nearest neighbor oxygen ions 
which form distorted octahedra. The average distances 
between the Er 3+ and O 2- ions in ErlOg and Er2C>6 



FIG. 1. (Color online) Top panel: Positions of the Er 3+ ions 
within SrE^CL, with the two crystallographically inequiva¬ 
lent Er 3+ sites shown in different colors. The blue box rep¬ 
resents a crystal unit cell of the Pnam space group. Bottom 
panel: a schematic of the magnetic interactions between the 
Er 3+ ions in one of the zig-zag chains. A four-particle cluster 
used to model the magnetization process is highlighted by the 
dashed box. 


octahedra in SrEr 2 C >4 are equal to 2.26 and 2.28 A, re¬ 
spectively, at 80 K. Crystal fields at the Rl and R2 sites 
have point symmetry C s . The Rl- and R2-sublattices are 
labeled below by the indices p from 1 to 4 {r p = r PiR i) 
and from 5 to 8 (r p = r p _ 4 , R2 ), respectively (see Fig. [TJ. 
Ions with the indices p = 1 and p = 3 (2 and 4, 5 and 
7, 6 and 8) are connected by the inversion operation and 
form a ladder with legs parallel to the c axis and shifted 
relative to one another by c/2. In general, a ladder can 
be considered as a zigzag chain where distances between 
the neighboring ions lying on the adjacent legs, d±, and 
on the same leg, d 2 , are slightly different, and d\ > d 2 . 
Each leg in a ladder has two neighboring legs belonging 
to other ladders, and the projection of the Rl and R2 
chains results in a honeycomb network (of edge-sharing 
distorted hexagons) in the ab plane, see Fig. [Q 

The SrY 2 04 :Er 3+ (0.5 at.%) sample used in the EPR, 
optical and magnetization studies was aligned using X- 
ray diffraction. It had the shape of a rectangular paral¬ 
lelepiped with the dimensions of 2.5 x 2 x 1 mm 3 with 
the edges along the crystallographic a, b and c axes, re¬ 
spectively. EPR spectra were measured with the X-band 
Bruker ESP300 spectrometer. The standard ER4102ST 
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rectangular cavity with the TE 102 mode was used. The 
studies were performed at 15 K, and the sample temper¬ 
ature was controlled using an Oxford Instruments ESR9 
liquid helium flow system. 

The magnetization measurements of the aligned 
SrY 2 04 :Er 3+ (0.5 at.%) single crystal were performed 
using a Quantum Design Physical Property Measure¬ 
ment System (PPMS-9) Vibrating Sample Magnetometer 
(VSM). 

Fluorescence of the crystals was excited with either a 
pulsed tunable dye laser (Littrow type oscillator and am¬ 
plifier, linewidth of about 1 A) pumped by the second or 
third harmonic of a Nd-YAG laser (LQ129, Solar LS) in 
the visible spectral range or a Ti:Sapphire tunable laser 
with linewidth of about 0.4 A (LX325, Solar LS) pumped 
by the second harmonic of a Nd-YAG laser (LQ829, So¬ 
lar LS) in the near-IR range. The spectra were analyzed 
with an MDR-23 monochromator. The fluorescence sig¬ 
nal was detected by a cooled photomultiplier (PMT-106 
or PMT-83) in the photon-counting regime. The studied 
SrY 204 :Er 3+ (0.5 at.%) crystal was kept in helium vapor 
at a temperature of 4.2 K. 

The polycrystalline sample of SrEr 2 04 used in the 
neutron scattering experiments was prepared from high- 
puritv starting materials SrCCU and Er 2 03 following 
Ref. |l|. The same sample was previously used for high- 
resolution powder neutron diffraction experiments [3j 
which have revealed the absence of any significant 
amount of impurities or chemical disorder in the sam¬ 
ples. The inelastic neutron scattering (INS) experiment 
was performed with the HET time-of-flight spectrometer 
at ISIS, STFC Rutherford Appleton Laboratory, United 
Kingdom in the temperature range 5 to 140 K. For the 
INS experiment the powdered samples were wrapped in a 
thin Al-foil and mounted inside a thin-walled cylindrical 
Al-can. The sample mount was cooled in a top-loading 
closed-cycle refrigerator with He-exchange gas and the 
INS data were collected using neutrons with incident en¬ 
ergies Ei between 24 and 120 meV, and for scattering 
angles between 3 and 135 degrees. The elastic resolution 
(full width half maximum) was 1.5, 2.0 and 5.5 meV for 
the Ei of 24, 40 and 120 meV, respectively. 

Although magnetization data for SrEr 2 04 have been 
previously reported SI3 the emphasis for the earlier 
studies was either on the magnetization behavior in high 
magnetic field [ 3 ! or on very low temperatures [l 6 j] where 
unusual plateaux have been observed rather than on the 
measurements of the absolute value of the magnetiza¬ 
tion. We therefore took advantage of the much higher 
accuracy available with a SQUID magnetometer (MPMS, 
Quantum Design) compared to the VSM and extraction 
magnetometers previously used [3] and remeasured mag¬ 
netization in SrEr 2 04 at T = 5 K for the field parallel to 
the three main crystallographical directions. The sample 
preparation is described in Ref. E3 



FIG. 2. EPR spectra (T = 15 K, v = 9.42265 GHz) of the 
SrY 204 :Er 3+ crystal for the magnetic field applied (1) along 
the c axis, (2) along the a axis, and (3) tilted by 20 degrees 
away from the a axis in the ab plane. 

III. EXPERIMENTAL RESULTS 
A. EPR spectra 

Examples of the observed EPR spectra taken in crossed 
static ( H ) and linearly polarized microwave ( H t ) mag¬ 
netic fields, in particular, for static fields applied along 
the crystallographic a and c axis, are shown in Fig. [2] 
For H || c and H t || b , two signals are detected in the 
spectrum at H ~ 680 Oe and ~ 1920 Oe. The first 
signal marked as Erl consists of the intense line originat¬ 
ing from erbium isotopes with even mass numbers and 
a pronounced hyperfine octet due to the single isotope 
with an odd mass number, 167 Er, which has a nuclear 
spin I = 7/2 and the natural abundance of 22.9%. The 
second signal at 1920 Oe marked as Er2 is considerably 
broader than the first one. As a consequence, it reveals a 
significantly less pronounced but still detectable hyper- 
fine structure of 167 Er 3+ . It should be emphasised that 
we cannot a priori identify the observed EPR signals with 
the R1 and R2 positions of the Er 3+ ions. The identifi¬ 
cation of the EPR signals (as well as of the transitions 
in the optical spectra discussed below) was achieved on 
the basis of the crystal field calculations (see Sec. m- 
The integral intensity of the resonance absorption is pro¬ 
portional to the square of the g-tensor component corre¬ 
sponding to the microwave field direction. The relation¬ 
ships between the intensities of the Erl and Er2 EPR 
signals (which are proportional to field derivatives of the 
absorption) in Fig. [2] agree with the ratios of the corre¬ 
sponding g-factors of the erbium ions at the R1 and R2 
sites determined below and presented in Table Q] 

With the magnetic field H applied along any of the 
crystallographic axes as well as lying in either the ac or 
be plane, two signals are observed in the EPR spectra in 
agreement with the presence of the two crystallographi- 
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TABLE I. g-factors of the ground doublets of the impurity 
Er 3+ ions in SrY 2 C >4 (the orbital reduction factors used in 
the calculations are k = 0.985 for Erl and k = 0.97 for Er2). 



Erl 

Er2 


Measured 

Calculated 

Measured 

Calculated 

5i 

5.382 ± 0.011 

5.396 

13.2 ± 0.1 

13.22 

92 

3.009 ± 0.007 

3.009 

1.69 ± 0.05 

1.726 

93 = <7cc 

9.93 ± 0.03 

9.944 

3.42 ± 0.02 

3.51 

<t> 

(48.0 ±0.2)° 

45.3° 

(8.1 ±0.05)° 

<p 

CO 

O 

Qaa 

4.239 

4.357 

13.07 

13.05 

9 bb 

4.478 

4.38 

2.502 

2.73 


cally inequivalent Rl and R2 sites. In an arbitrary ori¬ 
ented magnetic field tilted away from the ac (or be) plane, 
two magnetically inequivalent centers are revealed in the 
spectra for both the Erl and Er2 positions (see Fig. [2j 
spectrum number 3). 

In a low-symmetry crystal field, all the energy levels 
of an Er 3+ ion with the ground electronic configuration 
4 f 11 are Kramers doublets. In an external magnetic field, 
the effective spin-Hamiltonian (S = 1/2) of a Kramers 
doublet in the Cartesian system of coordinates with the 
coordinate axes along the principal axes of a g-tensor has 
the following form 

Us = + g2H 2 S2 + g^H^Sz). (1) 

Here mb is the Bohr magneton, H a and S a (a = 1, 2, 3) 
are projections of the magnetic field and effective elec¬ 
tronic spin moment on the coordinate axes, respectively. 
In the case of a local C s symmetry of a paramagnetic ion, 
two principal axes (labeled by the indices 1 and 2) of the 
< 7 -tensor are in the ab (mirror) plane, and the third axis 
is parallel to the c axis. 

In the EPR spectra with the frequency v and magnetic 
field lying in the ab plane at the angle $ to the a axis, 
the resonant fields H of erbium isotopes with even mass 
numbers are given by the equation 

, ( 2 ) 

) cos(2i? ± 2 ip K ) 

where +<p K and — tp K are the angles between the principal 
axes corresponding to the ^-factors gi tK and the a axis 
for the Er 3+ ions at sites r i jK , and r 2 , K , r^ K , respec¬ 

tively (k = Rl, R2). In agreement with the results of the 
measurements for an arbitrary orientation of the field H , 
four EPR signals can be observed, but when the magnetic 
field is applied along a principal crystallographic axis a 
(a = a,b or c), there are only two lines corresponding to 
the resonant fields H = 2nhu/(g aaiK iJ,B), where 

gaa, K = [Si,« + 52,« + (5i,« - 92, J cos(2<^ K )] 3 / 2 / v^, (3) 
9bb, K = \gl, K + gl, K - (5i,« - 92, K ) cos(2vj k )] 1/2 /v / 2, (4) 
9CC,K = 53 , K- (5) 


Mb H = 


2v27rhis 


\/ g\, K + &, K + {gl, K - 9l, K 


Hlla Hllb Hlla Hllc Hlla 



FIG. 3. (Color online) Measured (symbols) and calculated 
(lines) angular dependencies of the resonant magnetic field 
in the (a) ab and (b) ac planes of SrY 2 C> 4 :Er 3+ (0.5 at.%) 
at the frequency of 9.42265 GHz. Black/red solid symbols 
and gray/orange open symbols are used for the Erl and Er2 
centers respectively. 


For the magnetic field in the ac plane, the resonant fields 
satisfy the condition 

Mb H = 2t Thv[g 2 cc ^ + ( 5 2 q k - g 2 CCK ) cos 2 i?] _1/2 , (6) 

where $ is the angle between the field and the a axis. 

Angular dependencies of the EPR spectra were mea¬ 
sured for the magnetic field lying in the ab and ac planes. 
The resonant field dependence of the field direction for 
the erbium isotopes with even mass numbers are shown 
in Fig. [3] These dependencies are well described by equa¬ 
tions © and ©. The results of the fits to the experi¬ 
mental data, within the model of the isolated effective 
spin-1/2 doublet with the anisotropic g-tensor, are rep¬ 
resented in Fig. [3] by the solid lines. The values of the 
principal g-tensor components, g-factors in the directions 
of the crystallographic axes as well as the angles between 
the principal directions of the g-tensor and the crystallo¬ 
graphic axes are presented in Table Q] It should be noted 
that a possible misalignment of the sample in an external 
magnetic field may introduce errors in the obtained val¬ 
ues of the ^-factors of up to 0.1. Nevertheless, the EPR 
studies reported here provide conclusive evidence for the 
strong magnetic anisotropy of the Er 3+ ions of an easy- 
axis type with the easy-axis along the c direction at the 
Rl sites and along the two non-collinear directions in the 
ab plane at the magnetically inequivalent R2 sites. 

The available experimental data do not allow for a 
full characterization of the hyperfine interaction for any 
167 Er 3+ center. The values of the hyperfine constants 
can be given with a satisfactory accuracy for only a lim¬ 
ited range of the magnetic field directions with respect to 
the crystal axes. The restrictions arise due to the strong 
broadening of the EPR lines. 

Additional weak EPR signals were observed in the 
magnetic fields lying in the ac plane. The nature of the 
corresponding paramagnetic centers is not entirely clear, 
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FIG. 4. (Color online) (a) Luminescence spectra at T = 4.2 K 
of the R1 centers corresponding to the emission from the low¬ 
est sublevels of 4 S3/2 (lower spectrum) and 4 In/2 (upper spec¬ 
trum) multiplets into the sublevels of the ground 4 Iis /2 mul- 
tiplet induced by radiation with the wave numbers of 19193 
and 10355.5 cm -1 absorbed due to the excitations of the Er 3+ 
ions into the sublevel number 4 in the 2 Hn/ 2 multiplet or the 
sublevel number 6 in the 4 In /2 multiplet, respectively (see Ta- 
blc[TT]). (b) Luminescence spectra of the R2 centers induced by 
radiation with the wave numbers of 19273.6 and 10393.9 cm -1 
absorbed due to the excitations into the sublevels number 6 
in the 2 Hn / 2 and 4 In / 2 multiplets, respectively. The insets 
give graphical representations of the crystal field energy levels 
for the two crystallographic sites. 

but it is possible that the impurity Er 3 +ions substitute 
not only for Y 3+ but also for the Sr 2+ ions. 


B. Optical site-selective excitation and emission 
spectra 

We studied the emission of the SrY 2 04 :Er 3+ samples 
for two channels corresponding to the radiative transi¬ 
tions from the two metastablc states of the Er 3+ ions, 
namely, from the lowest CF sublevels of the 4 In /2 and 



FIG. 5. (Color online) Excitation spectra of the R1 (solid 
line) and R2 (dashed line) centers in SrY 2 04 :Er 3+ (0.5 at.%) 
corresponding to the transitions from the ground state into 
the 4 F 7 / 2 multiplet of the Er 3+ ions at T = 4.2 K registered 
by the emission from the lowest sublevels of the 4 In /2 mul¬ 
tiplet to the first excited sublevels of the ground multiplet 
at the wave numbers of 10170.7 cm -1 (10195.9 cm -1 - 25.2 
cm -1 ) and 10150 cm -1 (10191.1 cm -1 - 41.1 cm” 1 ), respec¬ 
tively. The arrows in the inset show the excitations (Exc), 
the registered emission (REm) and the nonradiative stepped 
decay of the excited states (NR.D). 


4 S 3/2 multiplets, to the sublevels of lower-lying multi¬ 
plets. 


The site-selective technique revealed the presence of 
two different sets of photoluminescent spectral lines orig¬ 
inating from transitions in the Er 3+ ions (see Fig. [4j. 
Here, energies of the CF sublevels of the ground ( 4 Ii 5 / 2 ) 
and the first excited ( 4 Ii 3 / 2 ) multiplets were determined 
from the selectively excited luminescence spectra. The 
selective excitation of different centres is confirmed, in 
particular, by the very different luminescence spectra 
shown in Figs. [Ur and (3 )d. 


Additional spectral lines were observed with an inten¬ 
sity that was at least two orders of magnitude weaker. 
These lines probably correspond to the Er 3+ ions substi¬ 
tuted for the Sr 2+ ions in the SrY 2 04 lattice, however, 
the concentration of these centers is very small. Ener¬ 
gies of all CF sublevels of the 4 In/ 2 , 4 Iq/ 2 , 4 Fg/ 2 , 2 H n / 2 , 
4 F7/2, 4 F5/2, 4 F 3 / 2 and 4 S 3 / 2 multiplets of both the R1 
and R2 optical centers were determined from the excita¬ 
tion spectra of the luminescence measurements (an ex¬ 
ample is shown in Fig. [5J. The measured spectra provide 
evidence of a relatively weak exchange of the excitation 
energy between the R1 and R2 centers. 
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TABLE II. Measured and calculated crystal field energies 
(in cm" 1 ) of the Er 3+ ions in SrY 2 0 4 :Er 3+ (0.5 at.%) and 
SrEr 2 0 4 (in brackets). 


Multiplet 

Erl 

Er2 


Measured 

Calculated 

Measured 

Calculated 

4 Il5/2 1 

0 

0 

0 

0 

2 

25.2 (26.2) 

25.25 

41.1 (41.1) 

42.6 

3 

75.1 

75.3 

74.6 (74.6) 

68.3 

4 

98.8 (97.3) 

97.3 

92.4 (93.4) 

94.2 

5 

111.2 (111.7) 

109.2 

162.2 (162.2) 

165.4 

6 

463.5 

457 

431.8 

427.6 

7 

491.6 (491.6) 

486.1 

453.6 (453.6) 

457.3 

8 

506 

497.4 

533.6 

528.4 

4 1-13/2 1 

6503.8 

6504.3 

6517.8 

6517.7 

2 

6540.5 

6540.2 

6568.1 

6560.6 

3 

6556 

6556 

6579.3 

6575.8 

4 

6580.5 

6576.3 

6635.1 

6632.3 

5 

6826.3 

6829.4 

6823.9 

6815.9 

6 

6842.1 

6838.4 

6831.6 

6822.9 

7 

6850.8 

6844 

6917.9 

6901.1 

4 Ill/2 1 

10196 

10197 

10192 

10193 

2 

10219 

10219 

10231 

10230 

3 

10232 

10231 

10258 

10259 

4 

10342 

10349 

10336 

10343 

5 

10348 

10354 

10340 

10346 

6 

10355 

10359 

10394 

10405 

4 l9/2 1 

12331.7 

12332 

12343 

12345 

2 

12404.9 

12400 

12430 

12425 

3 

12524.6 

12531 

12533 

12539 

4 

12581.4 

12584 

12584 

12599 

5 

12621 

12621 

12622 

12630 

4 F 9 /2 1 

15129.1 

15132 

15138.6 

15142 

2 

15156.1 

15156 

15176.6 

15186 

3 

15301.4 

15308 

15276.7 

15284 

4 

15333.8 

15340 

15377.9 

15386 

5 

15446.8 

15440 

15457.4 

15457 

4 s 3/2 1 

18184.8 

18190 

18201 

18205 

2 

18234.8 

18235 

18295 

18295 

2 H„/2 1 

19048.1 

19052 

19054.1 

19057 

2 

19066.3 

19063 

19072.1 

19079 

3 

19074.5 

19076 

19088.5 

19104 

4 

19193 

19207 

19205 

19213 

5 

19209.5 

19228 

19229 

19242 

6 

19265.2 

19271 

19273.7 

19274 

4 F 7/2 1 

20379 

20388 

20358.8 

20351 

2 

20437.5 

20438 

20462.2 

20452 

3 

20478.8 

20484 

20493.8 

20493 

4 

20543.9 

20535 

20585.8 

20575 

4 F 5 /2 1 

22047.2 

22063 

22053.5 

22061 

2 

22070 

22078 

22087.5 

22100 

3 

22194.7 

22179 

22193.8 

22182 

4 F 3 /2 1 

22451 

22453 

22428.8 

22428 

2 

22514.1 

22514 

22570.5 

22566 


Crystal field energies for both the R1 and R2 centers 
obtained from the analysis of the experimental data are 
presented in Table HT1 The obtained energy level patterns, 
and the total splittings of the ground and excited mul- 
tiplets are similar to those of the impurity Er 3+ ions at 
the Y 3+ sites with C 2 symmetry in Y 2 O 3 [Tsl j . The en¬ 
ergy level pattern of the Er 3+ ground multiplet 4 Ii 5/ 2 in 
the octahedral CF contains two groups of sublevels, sepa¬ 
rated by the energy gap of about 300 cm -1 in the case of 
oxygen ligands, the lower one involves a quadruplet Ts, 
a doublet r 6 and a second quadruplet T 8 , and the upper 
one involves a doublet T 7 and a quadruplet Tg, the total 
spitting being close to 500 cm -1 . Due to a splitting of 
the quadruplets by the low-symmetry component of the 
crystal field, one may expect to observe a pattern formed 
by the lower five doublets well separated from the upper 
three doublets. The measured CF energies for both the 
Erl and Er2 centers agree well with these simple argu¬ 
ments. 

The two optical centers formed by the Er 3+ ions in 
SrY 2 C >4 differ markedly in the splittings of the 4 S 3 / 2 
and 4 F 3 / 2 multiplets (see Table ITT1) . Large splittings of 
the multiplets with the total angular moment J = 3/2 
are likely be related to the more pronounced distortion 
of the first coordination shell (oxygen octahedron) for 
the Er 3+ ions at the R2 sites and, consequently, stronger 
quadrupole components of the crystal field. This con¬ 
jecture is confirmed below by the CF calculations (see 
Sec. HVAl) . 


C. Inelastic neutron scattering 

The results of the inelastic neutron scattering mea¬ 
surements on a powder sample of SrEr 2 04 are sum¬ 
marised in Fig. [Gj At a base temperature of 5 K, sev¬ 
eral peaks are clearly visible in the energy transfer range 
of up to 15 meV, they are labeled as A, B, C, D, and 
E. The widths of the peaks are resolution limited, apart 
from peak D, which consists of two excitation levels (see 
Sec. IIVBI) . The peaks show very little dependence (a 
slight decrease) of intensity with the scattering vector Q , 
consistent with the Q-dependence of the magnetic form 
factor for the Er 3 + ions. There were no detectable phonon 
signals observed at any Q, therefore the intensity of the 
peaks is presumed to be fully magnetic in nature. 

With increasing temperature, the main peaks start to 
look broader and several weaker peaks become visible 
at different energy transfers. This behavior is typically 
associated with the increase of the fractional population 
of the higher-energy levels at higher temperatures, which 
allows excited state transitions to become visible. 

Additional measurements with a higher incident neu¬ 
tron energy revealed the presence of a much weaker peak 
at about 20 meV (see Fig. 0) as well as a weak double 
peak at about 60 meV (not shown). 
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FIG. 6. (Color online) Neutron scattering spectra of a poly¬ 
crystalline SrEr 2 C >4 sample for the incident neutron energy 
Ei = 24 meV measured at different temperatures. The arrow 
points to the additional peak in the scattering intensity ob¬ 
served with the incident energy of 40 meV. The origin of the 
peaks labeled from A to E is discussed in the text. The results 
of simulations are represented by solid lines. The curves at 
different temperatures are offset for clarity. 



FIG. 8. (Color online) Temperature dependence of dc mag¬ 
netic susceptibility of SrEr 2 04 crystal measured along the 
three crystallographic axes (symbols) in a field of H = 1 kOe 
in the temperature range of 0.5 - 50 K [T^, [T?]. Solid lines 
represent the calculated susceptibility for temperatures from 
1 to 50 K in the framework of the four-particle self-consistent 
model (see Sec. II V Cl for further details). Inset: the tempera¬ 
ture dependence of the inverse susceptibility. 


D. Magnetometry studies 

Magnetization results at T = 2 K for a diluted 
SrY 2 04 :Er 3+ (0.5 at.%) sample are shown in Fig. [7] while 
Figures [5] and ED summarize the results of the suscepti¬ 
bility and magnetization measurements in the SrEr 2 04 
single crystals. 

For the diluted sample, the magnetization increases in 
an applied field in a way which preserves the order of the 
absolute values of magnetization observed in a concen¬ 
trated sample at higher temperatures, Mjj\\ a > Mh\\ c > 



FIG. 7. (Color online) Measured (symbols) and calculated 
(lines) magnetic field dependencies of the magnetization along 
the a, b and c crystallographic directions in the SrY 204 :Er 3+ 
(0.5 at.%) sample at T = 2.0 K. 


Mn\\b- The gradients of all three curves gradually de¬ 
crease with an applied held, as would be expected for a 
system of noninteracting magnetic moments. The details 
of the magnetization calculations for a diluted sample are 
given in Sec. IIV Al 

In Fig. [8] we combine the previously reported results 
for the temperature dependence of the susceptibility of 
SrEr 2 C >4 [16], U3| with the results of the calculations. For 
all three directions of an applied magnetic held, the four- 
particle model described in detail in Sec. IIV Cl returned a 
satisfactory agreement with the experimental data, par¬ 
ticularly at higher temperatures. A linear temperature 
dependence of the inverse susceptibilities along the three 
crystallographic axes is observed at T > 15 K (see in¬ 
set in Fig. [8]). The negative values of the corresponding 
Curie-Weiss temperatures [l6| indicate antiferromagnetic 
interactions between the Er 3+ ions. 

Fig. [9] shows a comparison between the experimen¬ 
tal and the theoretical results for the magnetization of 
SrEr 2 C >4 at T = 5 K. Although the M(H) curves have 
been reported previously, in particular in high helds Q 
and at very low temperatures [I(|, comparison with the 
more accurate experimental data helps to justify the va¬ 
lidity of the theoretical approach used (see Sec. IIV Cl for 
further details). One important observation to make here 
is that the magnetization for H || a is largely due to the 
Er 3+ ions at the R2 sites, while for H || b and H || c, 
the dominant contribution is from the Er 3+ ions at the 
Rl sites. Combined with the results of polarized neutron 
scattering [I|, which suggested that at low temperatures 
the magnetic moments in the long-range and short-range 
structures of SrEr 2 04 are predominantly pointing along 
the c and a axes, respectively, this observation uniquely 





























is the free-ion Hamiltonian written in a standard form 
that includes the energies of the spin-orbit interaction, 
electrostatic interactions between the 4/ electrons (la¬ 
beled by the index j), and additional terms due to inter¬ 
configuration and relativistic interactions, lj and §j are 
the orbital and spin moments of the electrons, L is the to¬ 
tal orbital moment, explicit expressions for the operators 
G, /, p, i. rh are given in literature (see Ref. 0 )- In the 
present work, we use parameters of the Hamiltonian (|5]) 
from Ref. [20j for the Er 3+ impurity ions in BaY 2 F§ mod¬ 
ified slightly to fit the gaps between the multiplet centers 
of gravity. Namely, the same parameters F 2 = 96320, 
F 4 = 67854, F 6 = 54260, a = 18.2, /3 = -600, 7 = 1760, 
P 2 = 594, P 4 = 445, P 6 = 119, M° = 3.7, M 2 = 2.07, 
M 4 = 1.15, T 2 = 390, T 3 = 50, T 4 = 83, T 6 = -271, 
T 7 = 298, T 8 = 280 (cm -1 ) for both the R1 and R2 sites, 
but different spin-orbit coupling constants ( = 2361 and 
2366 cm -1 are used for the R1 and R2 sites, respectively. 

The Hamiltonian Hcf gives the energy of the local¬ 
ized 4/ electrons in a static crystal field. In the crys¬ 
tallographic coordinate system with the z axis along the 
c axis, this Hamiltonian Hcf for a RE ion in the crystal 
field with the C s symmetry in SrP 204 , is determined by 
15 CF parameters B q : 

Hcf = B° 2 0° 2 + B 2 2 0 2 2 + B 2 2 0 2 2 + B°O 0 4 + B 2 0 2 A (9) 

+ B 4 2 0 4 2 + B 4 Oj + P 4 - 4 04 4 + B° 6 0° 6 + B 2 0 2 
+ B^ 2 Of + B 4 6 0 4 + Bq 4 Oq 4 + BlO% + P 6 - 6 0 6 - 6 . 


FIG. 9. (Color online) Measured (symbols) and calculated 
(solid lines) field dependence of the magnetization of SrEr 2 C >4 
single-crystals along the three crystallographic axes at T = 
5 K. The calculated magnetic moments at sites R1 and R2 are 
shown by dashed lines marked as Erl and Er2, respectively. 


identifies the R1 sites with the commensurate long-range 
order and the R2 sites with the incommensurate short- 
range order. 


IV. MODELING AND DISCUSSION 

A. Crystal fields affecting the Er 3+ ions at the R1 
and R2 sites 

The analysis of the experimental data presented in 
Sec. HU was carried out starting from the Hamiltonian, 
H, of a single Er 3+ ion operating in the total space of 
364 states of the ground electronic 4 / 11 configuration: 

H = Hfi + Hcf- (7) 

Here 

Hfi = C 'y ' ijSj + cxL “ + / 3 G{G 2 ) ^G{Ry) 
j 

+ Y / ( Fq f«+ pq Pi+ Tqi <i+ Mq ™'i) ( 8 ) 
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Here O are linear combinations of the spherical ten¬ 
sors Cq F [1I|, 0. In particular, O 3 = a p oCi P \ 0% = 
a pq (dj p) +C^), Op = -ia pq (Ci p) -C {p l) for q = 2,4,6; 
the numerical factors are 020 = 2, 022 = 2/x/6, Q 4 o = 8, 
042 = 4/VTO, 044 = 8/v/70 , 060 = 16 , 062 = 16/v/105, 
064 = 16/VT26, o 66 = 16//23I. There are two sets 
of independent CF parameters for the RE ions at the 
nonequivalent R1 and R2 positions, parameters B~ q (q > 
0) for the ions in the sublattices 1, 3 and 2, 4 (5, 7 and 
6 , 8) have the same absolute values but differ in sign. 

The initial values of the CF parameters were calculated 
in the framework of the semi-phenomenological Exchange 
Charge Model (ECM) 0 using the expressions 


B q = K q S^ — 
p p ^ R v 


Qi/(l o>) j^ p 


+ 2(2 p + !) s p (R v ) 


(r p J 

>P 

'V 

o*(e v ,$ v ), (io) 


where = o - Q 2 , I\ q = a~ 2 /2 for q / 0 , e is the ele¬ 
mentary charge, R u , 0^, are spherical coordinates of 
the ion located at the lattice site v, eQ „ are ion charges, 
cr p are the shielding factors, and ( r p ) are the moments of 
the 4/ electron density. Linear combinations S P (R V ) of 
the squared overlap integrals between the wave functions 
of the 4/ electrons localized at the RE ion and the elec¬ 
tronic wave functions of the lattice ions determine the 
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“exchange” charges. Because the overlap integrals de¬ 
crease exponentially with increasing distance between the 
ions, we take into account only the overlapping between 
the Er 3+ \4fl z ) wave functions and the |2s), |2 pl z ) wave 
functions corresponding to the outer electronic shells of 
the six nearest neighbor oxygen ions (here l z is the projec¬ 
tion of the electronic orbital moment on the bond direc¬ 
tion). In this case, introducing the dimensionless model 
parameters G s , G a and G n , we obtain [l3| 

S P (R„) = G s S s {R„) 2 + G a S a {R u ) 2 

+ ( 11 ) 

The R dependencies of the overlap integrals S s (R) = 
(4/0|2 a), S a (R) = (4/0|2p0) and S n (R) = (4/l|2pl), 
computed with the available radial wave functions of the 
Er 3+ [H and 0 2 ~ [H ions, can be approximated by the 
functions S u = a u exp(— b u R Cu ) (u = s,cr and 7 r), where 
a s = 0.333, b s = 1.0776, c s = 1.4269; a a = 0.08648, 
b a = 0.3845, Cfj = 1.9254; = 3.1067, b„ = 3.0381, 

c n = 0.8048 (R in angstroms). The results of the calcula¬ 
tions of the CF parameters (flOl) performed in the frame¬ 
work of the simplest single-parameter version of ECM 
(i G s = G a = Gn = 8 ) are presented in Table IIIII The 
values of (r 2 ) = 0.666, (r 4 ) = 1.126 and (r 6 ) = 3.978 
(in atomic units) (23[, og = 0.617 [25|, er 4 = a 6 = 0, 
the nominal values of the ion charges Q(Sr)=2, Q(Y)=3, 
Q(0)=-2 and the structure parameters of SrEr 20 4 [l5| 
were used in the calculations. It should be noted that 
the oxygen coordinates in the unit cell of SrY 20 4 deter¬ 
mined in Refs. EMU should be revised because they 
lead to significantly overestimated distances between the 
Y 3+ and the nearest neighbor O 2- ions as compared to 
the corresponding distances in SrEr 20 4 . In particular, 
according to Ref. the average Y 3 + -0 2_ distances 
in Y1C>6 and Y2C>6 octahedra exceed the corresponding 
average Er 3 +-0 2 distances in SrEr 2 0 4 by 0.08 A. 

The energy level schemes and the ratios between the 
3 -tensor components of the ground doublet for the Er 3+ 
ions at the R1 and R2 sites obtained from the numerical 
diagonalization of the Hamiltonian dTJ) which make use 
of the calculated CF parameters agree only qualitatively 
with the experimental data. In particular, the splittings 
of the 4 S 3 /2 and 4 F 3 / 2 multiplets are underestimated by 
about 30%. 

The next step in the simulations involves variations of 
the initial values of the CF parameters. To conserve the 
physical meaning of the CF parameters during the fit¬ 
ting procedure, namely, their definition within the crys¬ 
tallographic coordinate system, the CF parameters were 
varied to match simultaneously not only the measured 
energies of the CF levels (see Table HO) but the princi¬ 
pal values and orientations of the principal axes of the 
g-tensor of the ground doublet of the Er 3+ ion (see Ta¬ 
ble HJ) as well. 

When a strongly diluted paramagnetic crystal 
SrY 20 4 :Er 3+ is placed in an external magnetic field H , 
the Hamiltonian (Q of an Er 3+ ion contains an addi- 


TABLE III. Crystal field parameters Bp (cm -1 ) for the Er 3+ 
ions at ri.ru and r i,R 2 sites {Calc. - the results of calculations 
in the framework of the exchange charge model, Fin. - the 
final values from the fitting procedure). 



R1 

R2 

p 

g 

Calc. 

Fin. 

Calc. 

Fin. 

2 

0 

118.7 

188 

-33.1 

17 

2 

2 

118.6 

137.5 

-546.9 

-744 

2 

-2 

-104.2 

-171.2 

-110.8 

-125 

4 

0 

-59.8 

-57.3 

-55.9 

-60.2 

4 

2 

-1012 

-1066.2 

1064 

1033.2 

4 

-2 

1140 

1165.2 

-797.8 

-977.8 

4 

4 

-73.3 

-86.9 

379 

430.2 

4 

-4 

-1015 

-972.3 

-638.3 

-685.6 

6 

0 

-39.5 

-38 

-33.0 

-35.2 

6 

2 

-8.2 

-22.3 

-51.6 

-68.4 

6 

-2 

23.0 

22.8 

-13.6 

-42.8 

6 

4 

-8.6 

30.1 

-63.4 

-80.2 

6 

-4 

-145.6 

-115.2 

-186.8 

-191.4 

6 

6 

-208 

-162.2 

-66.7 

-119.6 

6 

-6 

-147.6 

-84 

56.3 

80.5 


tional term, the Zeeman energy, Rz = —rhH. Here 
rh = —ub + 2 §j) is the magnetic moment op¬ 

erator, the sum is taken over eleven 4/ electrons, and 
k is the orbital reduction factor that, for simplicity, is 
approximated by a scalar. The projection of the Zee- 
man energy on the 2 -dimensional space of wave functions 
(a = + and —) of a Kramers doublet T can be written as 
Rs — Mb Y2 HagapS/ 3 , where the effective spin S = 1/2, 

aft 

and components of the 3 -tensor are determined by the 
corresponding matrix elements of the magnetic moment: 

9 ax = 2 Re(T, +|m a |r, -), 

g ay = —2Im(T, +|m Q |r, —), 

9az — 2(r, +|m Q |r, +). (12) 

Square roots of eigenvalues of the matrix ( g 2 ) a /3 = 
X) 3 q 7 3 / 3 7 are the principal values of the 3 -tensor, and 

7 

eigenvectors of this matrix are the directional cosines of 
the corresponding principal axes. In the present work 
we considered the orbital reduction factor as an addi¬ 
tional fitting parameter with the condition that it may 
only deviate from unity by a small amount. The fit¬ 
ting procedure involved numerical diagonalization of the 
Hamiltonian 0 for each fixed set of the CF parame¬ 
ters and subsequent computations of the 3 -tensor com¬ 
ponents ( 1121 ) of the ground doublet; the main attention 
was paid to the correct description of the CF energies 
of the ground multiplet sublevels and principal values 
of the 3 -tensor. It is well established that the conven¬ 
tional one-electron CF parameterization underestimates 
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the splitting of the 2 Huj2 multiplet of Er 3+ [27l430| . Fol¬ 
lowing Refs. 0] and [30| . to fit the measured splittings, 
we introduced additional empirical parameters, the mul¬ 
tiplying factors of 1.37 and 1.32 for the Er 3+ ions at the 
R1 and R2 sites, respectively, in the reduced matrix el¬ 
ements of the fourth-rank spherical tensors in the space 
of states belonging to this multiplet. The final values of 
the CF parameters are presented in Table nm and the 
calculated ^-factors and the CF energies are compared 
with the experimental data in Tables Q] and ITT] respec¬ 
tively. Using the obtained sets of the CF parameters, we 
are also able to satisfactorily reproduce the results of the 
magnetization measurements. From a comparison of the 
calculated and measured field dependencies of the magne¬ 
tization along the crystallographic axes at a temperature 
of 2 K for the SrY 2 C> 4 :Er 3+ (0.5 at.%) sample assum¬ 
ing different populations c K of the Er 3+ ions at the R1 
and R2 sites (see Fig. 0. we obtained physically mean¬ 
ingful values of cri = 57% and cr 2 = 43%. This result 
proves the preferential occupation of the R1 sites by the 
impurity Er 3+ ions in the strongly diluted paramagnet 
SrY 2 0 4 :Er 3 +. 

B. Inelastic neutron scattering 

Intensity of neutron scattering from powder SrE^CU 
is considered here as a sum of two terms (see Eq. [13] be¬ 
low) corresponding to elastic and inelastic scattering. Ne¬ 
glecting magnetic interactions between the paramagnetic 
ions, i.e., supposing that the Er 3+ ions contribute inde¬ 
pendently to the scattering processes, we can write the 
scattering intensity versus energy transferred from the 
incident neutrons to the magnetic subsystem, av erag ed 
over directions of the scattering vector, as follows [3l| 

I(AE,T) = K f D(0,AE) (13) 

+ Ph,t{T)D{E k y' — E K r, AE)w k (T, T'), 

K rr' 

where 

D(x, A E) = (2tt 5 2 )~ 1/2 exp [—(a; - AE) 2 /2S 2 } (14) 

is the setup form-function that is approximated by a 
Gaussian, p Ki r(T) are populations of energy levels E K r 
with wave functions |acT, <t) of Er 3+ ions at K-sites, 
k =Erl, Er2; 

u> K (r,r')= ^ |(kT', cr'|TO a |rer, ct )| 2 (15) 

o’er', cx.=Xiy,z 

are relative probabilities of magnetic dipole transitions 
between energy levels of a RE ion, and Kf is a fitting 
parameter. Expression (ED is a reasonable approxima¬ 
tion for the scattering intensity assuming equal proba¬ 
bilities for the directions of the crystallographic axes in 
a powder sample. The spectral envelopes (1131) were cal¬ 
culated by making use of the transition energies given 


in Table El the transition probabilities computed with 
the wave functions obtained from the numerical diago- 
nalization of Hamiltonian 0 for impurity Er 3+ ions in 
SrY 2 04 , and a S v alue of 4.4 cm -1 . The corresponding 
HWHM of a/ 2 In 2 5 = 0.65 meV of the spectral lines is 
chosen to coincide with the average instrumental resolu¬ 
tion for the neutrons with an incident energy of 24 meV. 
The results are compared with the experimental data in 
Fig. [ 6 ] 

At low temperatures, in particular, at 5 K, the spec¬ 
trum is formed by transitions of the Er 3+ ions from the 
ground doublet to the excited crystal field sublevels of 
the 4 Ii 5/2 multiplet. According to the calculations, the 
peaks A and B (see Fig. © correspond to the transitions 
with integrated intensities (in units of the intensity of the 
transition A at the R1 sites) of 1.0 and 0.86 to the first 
excited doublets with energies 26.2 and 41.1 cm -1 at the 
R1 and R2 sites, respectively. Both the peaks C and D 
correspond to the unresolved closely spaced transitions to 
the third (with energies 75.1 and 74.6 cm -1 ) and fourth 
(with energies 97.3 and 93.4 cm -1 ) sublevels of the Er 3+ 
ions at the R1 and R2 sites, with the total intensities of 
0.45 + 0.30 = 0.75 and 0.35 + 0.26 = 0.61, respectively. 
The peak E with an intensity of 0.264 corresponds to 
transitions from the ground state to the fifth sublevel of 
the 4 Ii 5/2 multiplet with an energy of 111.7 cm -1 at the 
R1 sites. Note that the calculated intensities of tran¬ 
sitions from the ground state to the upper three sub- 
levels of the 4 Ii 5/2 multiplet at both the R1 and R2 sites 
are about two orders of magnitude weaker than the in¬ 
tensities of transitions within the lower groups of five 
sublevels. The temperature evolution of the spectrum 
envelope is caused by the redistribution of the sublevel 
populations, a number of additional transitions appear 
which lead to the smoothing of the scattering spectra at 
higher temperatures. 


C. Four-particle self-consistent model: dc magnetic 
susceptibility and the magnetization 

The above discussions prove that the Er 3+ ions at the 
R1 and R2 sites have rather different magnetic anisotropy 
in the ground state, the easy axis of the magnetization is 
parallel to the chain direction (along the c axis) for Erl, 
and is perpendicular to the chain direction for Er2. At 
temperatures of about 1 K the dipole-dipole interactions 
tend to induce a long-range ferromagnetic order along the 
chains at the Erl sites or an antiferromagnetic order of 
the Neel type at the Er2 sites because of small distances 
between the nearest neighbor ions in the chains and large 
values of the ^-factors <73 (Erl) and g\ (Er2). In particu¬ 
lar, when considering only the dipolar interactions, sim¬ 
ulations of the dc susceptibility along the c axis in the 
framework of the conventional single-site mean-field ap¬ 
proximation bring about a divergence at the Curie-Weiss 
temperature of 2.1 K. Powder neutron diffraction and 
single-crystal heat capacity measurements have demon- 



11 


strated long-range magnetic ordering in SrEr 204 below 
Tn = 0.75 K, but with a rather specific magnetic struc¬ 
ture: the magnetic moments at the Erl sites point along 
the c axis and form ferromagnetic chains along this direc¬ 
tion, with neighboring chains coupled antiferromagneti- 
cally [§:]. Note that the measured value of the magnetic 
moment at the Erl sites of 4.5/ig at T = 0.55 K is in 
line with the corresponding (/-factor g 3 = 9.93 which is 
determined in the present work. The measured magnetic 
moments of the Er 3+ ions at the R2 sites below Tn have 
much reduced values, and it has been shown that only 
short-range magnetic correlations emerge within the R2 
chains coexisting with the long-range order in the R1 
chains HI- 

In the following we focus our attention on the mag¬ 
netic properties of SrEr 204 in the paramagnetic phase. 
Each Er 3+ ion (at both the R1 and R2 sites) has eight 
neighboring Er 3+ ions at distances less than 0.4 nm. For 
example, an ion in the sublattice number 1 (at the R1 
site) has the two nearest neighbor ions belonging to the 
same sublattice (to the left and right within the chain, 
see Fig. [I]), a pair of the next nearest neighbor R1 ions 
belonging to the magnetically equivalent sublattice num¬ 
ber 3 (on the second leg of the ladder), and two pairs of 
more distant ions at the R2 sites belonging to the sublat¬ 
tices with numbers 6 and 7. There are good reasons to 
believe that the neighboring ions are connected by strong 
antiferromagnetic exchange interactions, such as the ob¬ 
served antiferromagnetic ordering of the R1 chains, and 
the measured temperature dependence of the suscepti¬ 
bility (see Fig. IS}, which demonstrates the suppressed 
responses of the Er 3+ ions in external magnetic fields at 
low temperatures. 

To elucidate the temperature dependence of the dc sus¬ 
ceptibility (in particular, the broad peak observed in the 
low-temperature susceptibility along the c axis), and the 
field dependence of the isothermal magnetization, we em¬ 
ploy the well known Bethe-Peierls approximation and in¬ 
troduce the four-particle self-consistent model. In the 
framework of this model we treat the dipole-dipole in¬ 
teractions exactly by computing the lattice sums with 
the Ewald method and consider the exchange coupling 
constants as fitting parameters. 

Let us consider four neighboring Er 3+ ions in a zig¬ 
zag chain labeled sequentially by indices I, II, III and IV 
(as shown in Fig. [[]). The effective Hamiltonian of this 
cluster has the following form: 

He =Uco-{mi+miv)Hi oc i-{rhii+rhiii)Hi oc2 , (16) 

where Hi oc 1 and if / oc 2 are local magnetic fields affecting 
the external and internal pairs of the ions (these two 
pairs of the ions - the first and the fourth, the second 
and the third - are not equivalent in the construction). 
The Hamiltonian Hco involves the single-ion energies Hj 
and the energies of interactions between the ions in the 


cluster: 

Hco = y ^ Hj + rhiJ r rhn + ran J'ram 

l=i, iv 

+ miuJrihiv + rhiJ c rhm + ranJ c mi V - (17) 

The matrices J r = Jr dd ^ + jl ex \ J' r = Jr dd ' > + 
and J c = J { c dd) + J determine the dipole-dipole (dd) 
and the exchange (ex) interactions between the Er 3+ 
ions along the rungs, r, and legs, c, of the ladder under 
consideration. The operators (fl6l) and cm work in the 
Kronecker product of the four single-ion Hilbert spaces 
of states. For magnetic fields of less than 60 kOe and 
temperatures below 60 K, computations of the magne¬ 
tization and susceptibility of the single Er 3+ ion at the 
R1 (or R2) site have shown that results obtained with 
the Hamiltonian determined in the total space of 364 
states of the electronic 4/ 11 configuration, and those with 
the single-ion Hamiltonian determined in the truncated 
Hilbert space spanned by the wave functions of the three 
lowest CF doublets, do not differ by more than a few 
percent. Based on these results, we constructed the ma¬ 
trix representation of the Hamiltonian (1161) in the basis 
of 6 4 = 1296 states corresponding to the three lower sub- 
levels of the ground multiplet of each of the four Er 3+ 
ions in the cluster. The contributions to the magnetic 
moments due to other excited single-ion states were taken 
into account in the last step of the simulations by adding 
the corresponding differences estimated from the exact 
single-ion calculations. 

To model the temperature dependence of the dc sus¬ 
ceptibility, first we compute the responses of the Er 3+ 
ions in the cluster to weak local magnetic fields. For 
a given set of the exchange coupling constants in the 
Hamiltonian m, components of single-site susceptibil- 

(i) 

ity tensors, X K ,ji are calculated numerically in the tem¬ 
perature range 1 — 50 K with a step size of 0.25 K. The av¬ 
erage magnetic moments of the Er 3+ ions are determined 
according to the following expressions (here, similar to 
the notation introduced above, k = Rl, R2; j = I — IV, 
(m Kj i) = (m K ,i V ), {m KtU ) = (m K , m )): 

{m Kjj ) = XkjHkJocI + XkJ H *,1 OC 2- (18) 

Such a cluster approach allows us to (at least partly) 
account for quantum correlations between the magnetic 
moments. Neglecting fluctuations of magnetic moments 
of the ions outside a fixed cluster (specifically, we consider 
a cluster containing the Er 3+ ions at magnetically equiv¬ 
alent Rl sites, belonging to the sublattices with numbers 
p = 1,3), we can write the local magnetic fields as fol¬ 
lows: 

H RUoc1 = ff R1 ,0 - 4 ex \m) m + AH m , (19) 

H m ,loc2 — H Rh0 + 4 dd \m) 

Rlj (20) 
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TABLE IV. Dipolar lattice sums Qpa,p '/3 (in units of 4n/3v c ). 


pa, p'/3 

Q P OL,p'p 

pa, p'l3 

Q poi,p' p 

px, px 

-4.116 

px, p + 4x 

-1.139 

py, py 

-4.781 

py, P + 4y 

4.144 

pz, pz 

11.897 

pz, 7z 

0.215 

lx, 2x 

0.629 

lx, 7y 

2.809 

1 y, 2 y 

2.370 

lx, 8x 

3.092 

lx, 3x 

-1.432 

i y, % 

-0.084 

1 y, 3 y 

3.717 

1 z, 8 z 

-0.008 

1 z, 3 z 

0.715 

lx, 8 y 

0.006 

lx, 3 y 

-4.657 

5x, 6x 

0.634 

lx, 4x 

2.310 

by, 6 y 

2.366 

ly, 4 y 

0.694 

5a:, 7x 

-1.471 

la;, 6x 

6.777 

5 y, 7y 

3.797 

i y, §y 

-4.082 

5 z, 7z 

0.674 

1 z, 6 z 

0.305 

5x, 7y 

-4.526 

1 ®, 6 y 

-0.045 

5a:, 8a; 

2.352 

lx, 7x 

-1.272 

5 y, 8 y 

0.653 

ly, 7 y 

4.058 

5 z, 8z 

-0.005 


where 

8 

H r1i0 = H~H d + Y^iQhP - J^Xrrip) 

p=i 

+ ( 4 dd) + 4 dd) - 4 ex) )(m) m . (21) 

Here H^, is the demagnetizing field, and (m) K are mag¬ 
netic moments of the Er 3+ ions at the sites k induced by 
the external field H (a more rigorous definition is given 
below). When H || c, only the z components of the 
average single-ion magnetic moments are nonzero, and 
{m PtZ ) = (m) R i, z for p from 1 to 4, (m PtZ ) = (m) R2 , z for 
p from 5 to 8. For the magnetic field H || a and H || b , 
the calculations of the local fields are more complicated 
because the single-ion susceptibility tensors have nonzero 
off-diagonal components, and the magnetic moments are 
tilted away from the external magnetic field direction in 
the ah plane. Specifically, for H || a, the components 
of the magnetic moments satisfy the relations (m PtX ) = 
(m) ri >x for p from 1 to 4, ( m PtX ) = {m) R2 ,x for p from 
5 to 8 , = -(m 2 , y ) = {m 3 , y ) = -(m 4j3/ ) = (m) m ,y 

and (m 5j y) = -(m 6 , w ) = (m 7 , y ) = -(m^ y ) = (m) R2 ,y 
Similar relations, with the rearranged indices x and y, 
are valid for if || b. The matrices J p e p ) in Eq. (1211) cor¬ 
respond to the exchange interactions between the next- 
nearest-neighbors and between the Erl and Er2 ions. 
The nonzero values of the dipolar lattice sums 


Q P 


E 

L(n 1 ,n 2 ,n 3 ) 


( r Lpp'^aP + 3Xl pp ' aXLpp' ()) 


( 22 ) 


Lpp' 


(here a'Lpp'a. is the a-component of the vector rL pp , = 
n±a + n 2 b + n 3 c + r p — r p >, and n\,n 2 ,n 3 are inte¬ 


gers) are presented in Table IIVI in units of 47r/3w c , 
where v c = abc is the unit cell volume). According 
to the crystal lattice symmetry Q 2 ,3 = Qi A ; <52,4 = 
<5l.3! <53,4 = <5l,2! <56,7 = <55,8! <56,8 = <5s,7; <57,8 = 
<55,6; <52,5 = <53,8 = <54,7 = <5l,6! <52,7 = <53,6 = 
<54,5 = Qi,s; <52,8 = <53,5 = <54,6 = <5 i, 7- The ad¬ 
ditional magnetic field Aif R i introduced in Eq. (1191) is 
determined from the condition of the magnetic equiva¬ 
lence of all ions in the cluster, namely 


(wi K , i) = {m K , h) = (m Ki in) = (m K) i V ). (23) 


This condition renormalizes the single-site susceptibili¬ 
ties, which now take the following form: 


v (») -v^fv^-yl 1 ) rV 2) 

X.K, — A.K,lfA. K ,I X K ,II) A.k,II 

- v (1) (v (1) - v (1) rV 125 
A. K .IllA. K ,I Xk,u) Xk.i- 


(24) 


The components of the bulk dc susceptibility tensor 
X = N a (x R i + AR2)/2 (per mole of Er, where N A is 
the Avogadro constant) are obtained from solutions of 
a system of self-consistent linear equations which deter¬ 
mine the susceptibilities Xk, of the Er 3+ ions at the Rl 
and R2 sites: 


(m) R i = xmH = XmHm,ioc2{{m) m , (m) R2 ), (25) 

{m) R2 = XR 2 H = X R2 -H'R2,ioc2((m) R i, (m) R2 ). (26) 


In the temperature range of 1—50 K, the computed sus¬ 
ceptibilities describing the response of SrEr 2 0 4 samples 
to external magnetic fields directed along the principal 
crystallographic axes agree well with the experimental 
data, as shown in Fig. [HJ Relatively small differences be¬ 
tween the measured and simulated susceptibilities along 
the a and b axes become noticeable only at temperatures 
below 3 K. It should be noted that the broad maximum 
observed in the susceptibility along the a axis (at tem¬ 
peratures close to 7 n = 0.75 K) is shifted in the calcu¬ 
lated temperature dependence of Xaa to 1.25 K. The re¬ 
sults of calculations hint that the low-temperature broad 
peaks in \aa and Xbb are related to the intrinsic properties 
of the Er2-ladders (he., the geometrically frustrated ex¬ 
change interactions), and that these peaks could perhaps 
be observed more clearly if the temperature range was 
extended to below 0.5 K despite the magnetic ordering 
of the Erl-ladders. The parameters of the exchange and 
dipolar interactions between the Er 3+ ions in the clusters 
used in the calculations are presented in Table |V] (note 
that in our model J r = J' r ). From the data, it can readily 
be seen that the introduced anisotropic exchange inter¬ 
actions are comparable with the dipolar ones. The addi¬ 
tional contributions to the local magnetic fields due to the 
exchange interactions (defined by Eq. I21T) are estimated 
from the fitting procedure. Namely, they are determined 

by the parameters £ = 13.4, £ j[ ex) = 

,9=1,3 p=l,3 

T(eaO _ j(ex) _ o q j(ex) _ 

Z-/ y,py J 5z,pz o.tf, J lz,pz 

p— 5,7 ,9=5,7 ,9=1,3 


14.5 
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TABLE V. Parameters of the exchange ( Jr ex \ ) and 

dipole-dipole ( Jr dd \ ji dd ' > ) interactions (in units of 1CP 3 
cm^ 1 /Mb) between the Er 3+ ions in the R1 and R2 ladders. 



R1 

R2 

j(-) 

J ot@ 

(ex) 

(dd) 

(ex) 

(dd) 

j(-) 


2.5 

4.1 

28 

4.2 

j(-) 

J r,yy 


2.5 

-7.1 

7.5 

-7.2 

j(-) 

>-'r,xy 

• a 

ll 

0 

10.3 

0 

10.0 



24.0 

3.0 

7.5 

3.0 

j(.) 

c,xx 

— Jc,yy 

5.0 

11.2 

5.0 

11.2 

j(-) 

c,zz 


5.0 

-22.4 

5.0 

-22.4 


r(ex) _ T (ex) _ 17 O T ( ex ) 

2-, x.px ~ L-* x.px — ^ J 5x.px 
p— 5,8 p—1,4 p= 5,7 

the units of 10 -3 cm _1 /^t|). 


31.2 (in 


At low temperatures, the energy of interaction deter¬ 
mined by the parameter is actually a bilinear 

function of the corresponding ^-factors of the ground 
state. Remembering the large value of g xx of the Er 3+ 
ions at the R2 sites, we conclude that the peculiar mag¬ 
netic properties of these ions (and particularly, the ab¬ 
sence of long-range order) are caused by strong antifer¬ 
romagnetic interactions between the ions forming edge¬ 
sharing triangles in the R2 ladders. In contrast, in the 
R1 ladders, a ferromagnetic dipolar interaction along the 
legs prevails due to the large value of the g -factor g zz , 
and induces the long-range magnetic order seen at low 
temperatures. 


The reliability of the parameters presented above and 
in Table [V] is confirmed by the good agreement between 
the calculated and measured field dependencies of the 
magnetization of SrEr 204 for the magnetic fields directed 
along different crystallographic axes as shown in Fig. [9] 
Calculations of the average self-consistent magnetic mo¬ 
ments at 5 K were performed using density matrices cor¬ 
responding to the Hamiltonians (fl6l) of the four-particle 
clusters in the R1 and R2 ladders and the procedure of 
subsequent approximations. The results of calculations 
(in Fig.© clearly demonstrate that there is a strong mag¬ 
netic anisotropy of the easy-axis type along the c axis at 
the R1 sites and along the a axis at the R2 sites. 


V. SUMMARY 

To summarize, we report the results of a system¬ 
atic study of the behavior of the magnetic Er 3+ ions 
in SrEr 2 04 as well as in a lightly doped nonmagnetic 
analogue, SrY 2 04 :Er 3+ . The observed EPR spectra in 
the single crystals of SrY 2 04 :Er 3+ (0.5 at.%) established 
very anisotropic ^-factors for the Er 3+ ions in the two 
crystallographically inequivalent sites. The CF energies 
of the Er 3+ ions substituting for the Y 3+ ions were de¬ 
termined by making use the methods of the site-selective 
laser spectroscopy. This combination of the EPR and 
spectral studies was further complemented by inelastic 
neutron scattering as well as the measurements of the 
field dependencies of magnetization in the paramagnetic 
phase of the concentrated system, SrEr 2 04 . The ob¬ 
tained data allowed us to develop a theoretical model 
to describe the electronic structure of the Er 3+ ions both 
in the magnetically dilute and concentrated samples. Al¬ 
though operating with a large number of independent 
parameters it was not possible to determine all of their 
values unambiguously, we are able to estimate most of 
these parameters from the literature data or from simu¬ 
lations based on known models developed for CF theory. 
The sets of CF parameters obtained for the Er 3+ ions 
at the two crystallographically inequivalent positions are 
related to the crystallographic axes, contrary to the re¬ 
sults of the CF simulations for the impurity Eu 3+ ions in 
Sri ? 2 04 (R = Y, Gd, In) reported in Ref. [3^ |. 

It is rather important to be able to separate the single¬ 
ion effects from the collective interactions for the over¬ 
all understanding of the complex low-temperature mag¬ 
netic properties of the concentrated strontium rare earth 
oxides SrR 2 04 (as well as the related family of barium 
rare earth oxides [33l - l35| b We hope that the approach 
adopted in this paper and that has proven successful for 
SrEr 2 C> 4 , could be applied to other compounds. In par¬ 
ticular, our initial simulations show an Ising type mag¬ 
netic anisotropy along the c axis and close to the b axis 
in the ground quasidoublet states of the Ho 3+ ions at 
the R1 and R2 sites, respectively, in SrHo 2 C> 4 . A more 
detailed analysis of the spectral and magnetic properties 
of the Ho 3+ ions in SrY 2 04 and SrHo 2 C >4 is in progress. 
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